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ABSTRACT 



Context. The interstellar medium contains warm clouds that are embedded in a hot dilute gas produced by supernovae. Because both 
gas phases are in contact, an interface forms where mass and energy are exchanged. Whether heat conduction leads to evaporation 
of these clouds or whether condensation dominates has been analytically derived. Both phases behave differently dynamically so that 
their relative motion has to be taken into account. 

Aims. Real clouds in static conditions that experience saturated heat conduction are stabilized against evaporation if self-gravity and 
cooling play a role. Here, we investigte to what extent heat conduction can hamper the dynamical disruption of clouds embedded in a 
streaming hot plasma. 

Methods. To examine the evolution of giant molecular clouds in the stream of a hot plasma we performed two-dimensional hydro- 
dynamical simulations that take full account of self-gravity, heating and cooling effects and heat conduction by electrons. We use the 
thermal conductivity of a fully ionized hydrogen plasma proposed by Spitzer and a saturated heat flux according to Cowie & McKee 
in regions where the mean free path of the electrons is large compared to the temperature scaleheight. 

Results. Significant structural and evolutionary differences occur between simulations with and without heat conduction. Dense 
clouds in pure dynamical models experience dynamical destruction by Kelvin-Helmholtz (KH) instability. In static models heat con- 
duction leads to evaporation of such clouds. Heat conduction acting on clouds in a gas stream smooths out steep temperature and 
density gradients at the edge of the cloud because the conduction timescale is shorter than the cooling timescale. This diminishes the 
velocity gradient between the streaming plasma and the cloud, so that the timescale for the onset of KH instabilities increases, and the 
surface of the cloud becomes less susceptible to KH instabilities. The stabilisation effect of heat conduction against KH instability is 
more pronounced for smaller and less massive clouds. As in the static case more realistic cloud conditions allow heat conduction to 
transfer hot material onto the cloud's surface and to mix the accreted gas deeper into the cloud. 

Conclusions. In contrast to pure dynamical models of clouds in a plasma and to analytical considerations of heat conduction that can 
evaporate such clouds embedded in a hot plasma, our realistic numerical simulations demonstrate that this destructive effect of KH 
instability is significantly slowed by heat conduction so that clouds can survive their passage through hot gas. 

Key words. ISM: clouds - ISM: structure - Conduction - Hydrodynamics - Instabilities - Method: numerics 



1. Introduction 



*h ' 1.1. The multi-phase Interstellar Medium 

The Interstellar Medium (ISM) is frequently described as an in- 
homogeneous ensemble of three phases (Mc Kee & Ostriker 
1977: IMQ77I) . The cold neutral phase with temperature T ~ 
80 K and density n ~ 40 cm -3 is represented by the cores 
of molecular clouds which are confined by a warm neutral to 
slightly ionized medium (T ~ 8000 K, n ~ 0.3 cm -3 ). These 
two components are in pressure equilibrium if the gas is ex- 
ternally heated and can cool radiative ly (Field et al. 1969). 
According to IMQ77I they are embedded in a third phase of the 
ISM: the hot dilute intercloud medium (HIM) with T ~ 10 6 K 
and n ~ 10~ 3 — 10~ 4 cm -3 . Since this component is produced 
locally by supernova explosions at even higher temperatures and 
densities, originally it cannot be in pressure equilibrium with the 
cooler phases and therefore has to expand. During this expansion 
shocks arise and the HIM penetrates the ambient clumpy ISM. 
Denser clouds cannot be swept-up by the shock front but are 



Send offprint requests to: G. Hensler, 

e-mail: hensler@ast.ro . univie . ac . at 



overrun and become engulfed by the HIM. Because of hot gas 
expansion the relative velocity between both gas phases leads to 
strong distortion of the clouds, resulting in a stripped-off gas tail. 
Clouds of this type can be found as remnants of larger clouds in 
galactic chimneys such as the one associated with the Hll re- 
gion W4 peyeret al. 1996| [Taylor etal. 1999| l or behind shock 
fronts of supernovae. Due to the strong discrepancy in the phys- 
ical states between the phases an interface has to form where the 
hot phase and the molecular cloud are in contact. Temperatures 
and densities are connected through steep gradients that lead to 
energy and mass transfer. 

This situation occurs in many astrophysical phenomena such 
as High-Velocity Clouds (HVCs), Hi structures characterized 
by radial velocities that are incompatible with simple models 
of the differential rotation of the galactic disk (see Wakker & 
van Woerden (1993 f° r a recent review). Interferometer mea- 
surements of the 21 cm Hi line at 1' resolution show several 
small clumps embedded in larger emission regions (Wakker & 
Schwarz [1991| l. This, along with their line width, led Wakker & 
Schwarz to conclude that the HVCs have a multi-phase struc- 
ture consisting of a cold, dense core and a warmer, more ten- 
uous halo. Distance measurements of the cloud complexes re- 
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main difficult. For at least two of them upper limits for their dis- 
tances are deduced by Danly et al. (119931 1, Keenan et al. ( 1995) 
and van Woerden et al. (1997 ) which place them in the galac- 
tic halo. From absorption line and X-ray observations it is well 
established that the galactic halo is filled with hot gas with a 
scaleheight of 4 kpc (IPietz etal. 19981 ) so that the HVCs have to 
interact with the hot gas. ROSAT observations show an increase 
of X-ray emission at the edges of the clouds that can be assumed 
as evidence for the existence of a hot interface (Kerp et al. 1 19941 
119951 ). The detection of velocity bridges ( IPietz et al. 19 96), con- 
nections of Hi gas in velocity-space-diagrams, could also stem 
from the interaction of HVCs with low or intermediate-velocity 
gas. Although the origin of most of the complexes is still specu- 
lative, HVCs are a classical example of cold multi-phase struc- 
tures moving through a hot plasma. 

Another possibility for the scenario described above ismuch 
more massive and larger proto-globular cluster clouds (PCCs) 
with temperatures near 10 4 K and densities several hundred 
times that of the surrounding gas. They therefore are gravita- 
tionally unstable at Jeans masses larger than 10 6 Mq. An upper 
limit for their mass can be approximated by the critical mass for 
an isothermal sphere that is embedded in a surrounding medium 
with pressure Pism (Ebert 1955; Bonner [19651 ): 

M max =i. 2 (M:) 2 G -^P I -v 2 (1) 

where G is the gravitational constant, ks Boltzmann's constant, 
T a PCC temperature and /i its mean molecular weight. These 
clouds are assumed to originate from condensations of thermally 
unstable gas with temperatures of some million Kelvin in the 
early epoch of galaxy formation and can be envisaged as pro- 
genitors of globular clusters (Fall & Rees 1985 ). The PCCs had 
to resist their gravitational collapse until star formation ignites 
with a very high efficiency. This condition is not easily fulfilled 
because the PCC is accelerated by the gravitational potential of 
the forming protogalaxy. While moving through the hot halo it 
will become subject to the growth of Kelvin-Helmholtz (KH) 
and Rayleigh-Taylor (RT) instabilities (Drazin & Reid 198Tb. 

1.2. Numerical models of clouds in a hot plasma stream 

To investigate an influence of KH and RT instabilities on the 
evolution of molecular clouds several authors have used numer- 
ical simulations. Murray et al. (1993) compared the evolution 
of clouds embedded in the subsonic stream of a dilute medium 
with and without self-gravity. Radiative losses are neglected and 
the gravitational potential of the clouds remain fixed to the ini- 
tial value of the calculation. Models without self-gravity and 
cloud masses M c \ <C M max break up after only a few dynam- 
ical times Tdyn = i? c l/ c cb where R c \ is the cloud radius and 
c c i is the sound speed inside the cloud. The mass loss after 2.5 
T(j yn is 20% and increases to 75% after 3.8 T^ yn . On the other 
hand, gravitationally bound clouds with a gravitational acceler- 
ation g = g cr i t , the value where perturbations with wavelengths 
of the order of R c [ are totally damped, show a different evolu- 
tion. The value for g cr [ t results from the relative velocity U, the 
cloud density p c \ and the density of the background material p^a 
to g crit = 2vr[/ 2 p bg /(p cl i? cl ). 

The gravity is sufficient to stabilize the cloud although the cloud 
edge is distorted. After 3.2 T<iy n only 2% of the initial cloud mass 
is lost and after 10 r^yn, 11%. Severing ([1995) improved the 
simulations done by Murray et al. ( 1993 ) by solving the Poisson 
equation for self-gravity self-consistently at each timestep and 



by including heating and cooling effects. The cooling takes a col- 
lisionally dominated plasma in thermal and ionization equilib- 
rium into account (Bohringer & Hensler 1989 ) and is balanced 
by heating for a cloud at rest. He demonstrated that clouds in a 
stream with Mach number 0.1 and with g > g CI [ t gain mass while 
at higher but still subsonic Mach numbers and with g < <ji cr it 
similar results to those by Murray et al. are achieved. Dinge 
(119971 ) applied cooling only to the gas ablated from the cloud 
and showed that cooling tends to accelerate the destruction of 
the cloud by RT instability. However, this arises due to the fact 
that the stripped cloudlets hit the cloud from the rear and thus 
trigger soundwaves which move through the cloud and stretch it 
along the symmetry axis. Globally, his results agree qualitatively 
with those of the former authors. 

Vietri et al. (1997 ) examined the influence of radiative losses 
on the evolution of KH instabilities. By analytic approximations 
they proved that plane parallel identical fluids are destabilized if 
they stream relative to each other at high Mach numbers while 
low Mach number flows tend to be more stable. For fluids having 
different densities, cooling processes exacerbate the KH instabil- 
ities and all Mach numbers become unstable although with mod- 
erate growth rates. For clouds without self-gravity this means 
that KH instabilities cannot be suppressed. The instability is con- 
strained to a small volume around the surface because cooling 
and heating timescales are shorter than the dynamical ones so 
that the instability cannot extend much beyond the interface. 
In addition, observationally derived structure parameters allow 
the conclusion that HVCs are bound by a significant Dark Matter 
halo mass. Quilis & Moore (2001) studied with hydrodynamical 
simulations the behavior of HVCs passing the dilute hot halo 
gas of a typical disk galaxy at large distances. The main purpose 
of their numerical models was to compare the head-tail struc- 
ture with observations to study the necessity for a Dark Matter 
content and a possible lower density limit of the surrounding in- 
tergalactic medium (IGM). Although they do not trace the evolu- 
tion with respect to mass loss and structure survival (e.g. the pure 
gas clouds are shown after 3 Myr only), the main conclusions 
are two-fold, namely, at first that the observed head-tail struc- 
tures with sufficient densities are produced if the wind density 
of the IGM exceeds iiicm > 10~ 4 cm~ 3 due to ram-pressure 
stripping (RPS) but independent of the existence of a DM halo. 
This supports the galactic fountain scenario for their origin be- 
cause the outermost IGM is expected to fall short of this density 
threshold. Secondly, the velocity was varied between 200 and 
400 km s _1 without any substantial change in the structure. 

1.3. Heat conduction 

Stable models for Giant Molecular Clouds (GMCs) like HVCs 
and PCCs assuming hydrostatic and thermal equilibrium consist 
of large temperature and density gradients at the surfaces of the 
clouds where the energy densities of the surrounding ISM and 
the cloud become equal. 

There are indications that the observed Hi masses of some HVCs 
are orders of magnitude too low to provide gravitational stabi- 
lization, so they must be confined by external pressure (Konz et 
al. |2002| ). On the other hand, the core-halo structure seen in some 
HVCs (e.g Wolfire et al. 1995) may result from minihalos com- 
posed of gas and dark matter that would likely move together as 
the clouds fall onto the Milky Way (Blitz et al. |1999l ). In this pic- 
ture, the largest HVCs are gravitationally bound by the dark mat- 
ter rather than pressure confined (Braun & Burt on 19991 2000). 
Because of the large temperature gradients and the high temper- 
ature of the surrounding medium of the order of some million 
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Kelvin, heat conduction must play a substantial role in the evo- 
lution of such clouds. 

Analytical studies of the influence of heat conduction on the 
evolution of clouds at rest were undertaken by several au- 
thors: Cowie & McKee (1977) (hereafter: ICM77I) investigated 
the evaporation of spherical clouds in general. For the classical 
regime described by a collision-dominated plasma they used the 
conductivity of Spitzer (1962). In this case, the heat flux due to 
heat conduction can be written as 



<?class = ~« • VT 

with 

1.84- 1CT 5 T 5 / 2 

K = 

In* 

where the Coulomb logarithm is 



erg s 1 K 1 cm 1 



In* = 29.7 + In 



T, 



6,e 



(2) 



(3) 



(4) 



with the electron density n e (in cm~ 3 ) and the electron temper- 
ature T 6 , e in units of 10 6 K. 

In a collisionless plasma the diffusion approximation for the heat 
flux breaks down because the mean free path of electrons be- 
comes comparable to or even larger than the temperature scale- 
height. In this so-called saturated regime described by a colli- 
sionless plasma, CM77 used a flux-limited heat flux. This takes 
charge conservation into account and yields results in good 
agreement with more sophisticated treatments (e.g. Max et al. 
1980) and with numerical simulations of laser heated plasmas 
(Morse & Nielsen [1973] Mannheimer & Klein |1975ll. T he satu- 
rated heat flux as an upper limit takes the form (ICM771 I 



|<Zsat| = 5<S> s pc J 



(5) 



with the sound speed c. $ s is an efficiency factor less than or of 
the order of unity, which embodies some uncertainties connected 
with the flux-limited treatment and flux suppression by magnetic 
fields (<f> s = 1 in our calculations). 

For an abrupt change of the conductivity from the classical to 
the saturated regime this leads to an envelope around the cloud 
consisting of three layers: a saturated zone embedded in an inner 
and outer classical zone. CM77 obtain a classical mass-loss rate 
of 



m class 



l6lTflKfR 

25k B 



(6) 



where Kf is the conductivity evaluated for the unperturbed hot 
medium far away from the cloud. R is the cloud radius and /i 
the mean molecular weight. The mass-loss rate is lower for the 
saturated case. In McKee & Cowie d 19771 ) they included radia- 
tive losses in their studies. As a criterion to separate the classical 
and the saturated case for a cloud of radius R embedded in a hot 
gas with temperature Tf and electron density n e /,|CM77] intro- 
duced a global saturation parameter q-q which is the ratio of the 
electron mean free path to the cloud radius 



0.08k/!/ _ 1-23 • 10 4 T| 



<f> s PfC 3 f R 



i ef 



R 



(7) 



For do < 0.027/$ s material condenses onto the cloud be- 
cause radiative losses exceed the conductive heat input. For 



0.027/$ s < cr 
tion, while for ao 
Begelman ( 1990) found similar results introducing the so-called 
Field length 



< 1 the cloud suffers classical evapora- 
> 1 the evaporation is saturated. McKee & 



Af = 



kT 



n 2 C 



M 



1/2 



with Cm = max(A,r/n) 



(8) 



in which the cooling or heating is comparable to the conductive 
energy exchange. For R > Ap condensation occurs, otherwise 
the cloud will be evaporated, n is the particle density, A the cool- 
ing rate, V the heating rate. 

These studies concluded that cold clouds in a dilute medium at 
a temperature of some million Kelvin always experience evap- 
oration, but none of them examined more realistic multi-phase 
clouds or included self-gravity. 

In a seperate paper (Vieser & Hensler 2005; hereafter: Paper l\ 
we studied the differences caused by heat conduction on the eva- 
poration/condensation competition between the fixed analytical 
description and the more realistic approach of flux saturation that 
adopts flexibly to the temporal physical state. 
The main results of these investigations are as follows: The an- 
alytical mass loss rates of a cloud at rest in a hot and rarefied 
medium can be reproduced in numerical simulations for the pure 
classical case, because the evaporated material is pushed away 
with supersonic speed. The initial large density and temperature 
jump at the edge of the cloud remains unaltered during the calcu- 
lation. Taking the more realistic saturated heat flux into account, 
a transition zone forms at the cloud edge in which the steep 
temperature and density gradients are reduced. This results in 
a lower evaporation rate than predicted. Simulations that include 
additional heating and cooling show an even more dramatic ef- 
fect. The clouds can even gain material if radiative cooling ex- 
eeds the energy input by heat conduction. 
Here we examine the evolution of molecular clouds in the stream 
of a hot, dilute medium. The treatment of heat conduction in 
the context of hydrodynamical simulations is described in §2. 
Analytical estimates of the influence of heat conduction are com- 
pared with the results of dynamical models with and without heat 
conduction in §3. Conclusions are drawn in §4. 

2. Hydrodynamical treatment 

The evolution of clouds in the subsonic stream of a hot plasma 
is studied by two-dimensional hydrodynamic simulations. 
The hydro-part of this Eulerian, explicit code is based on the 
prescription of Rozyczka (1985) and has been extensively tested 
and used by different authors (e.g. Yorke & Welz 1996). The hy- 
drodynamic equations have been formulated in cylindrical coor- 
dinates (r, z), assuming axial symmetry around the z-axis that is 
also the flow direction. The cloud's center is located on the z- 
axis. The differencing scheme used to discreticize the equations 
is second-order accurate in space because a "staggered grid" is 
used . We applied operator splitting for time integration, be- 
cause numerical experiments have shown that a multi-step so- 
lution procedure is more accurate than a single integration step 
based on preceding values flStone & Norma n"! 992| l. The advec- 
tion scheme of van Leer (1 19771 1 is employed. Since the basic 
code is explicit, the Courant-Friedrichs-Lewy (CFL) condition 
determines the maximum time step for the hydro-part. Because 
the conduction time step is smaller than the CFL one, the temper- 
ature distribution has to be calculated several times in one hydro 
time step. Von Neumann-Richtmyer artificial viscosity is used 
for the treatment of shocks. In order to prevent the cloud from 
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moving out of the computational domain due to drag forces, the 
cloud center of mass is re-adjusted at each time step. 
The grid parameters, the resulting physical domain and the res- 
olution are listed in Table [TJ for three representative models. The 
grid resolution is 28 - 33 zones per initial cloud radius. 



Table 1. Numerical parameters for the simulations 



numerical grid 



Model celr 



physical grid 
pc 2 



resolution 
pc cell" 1 



U 640 x 200 800 x 250 1.25 

K 640 x 200 80 x 25 0.125 

E 640 x 200 50 x 15.625 7.8 ■ 10' 



The boundary conditions on the upper and right-hand sides are 
semi-permeable to allow for an outflow of gas from the com- 
putational domain. The physical parameters at the lower bound- 
ary, the symmetry axis, are mirrored. The parameters, density 
and temperature at the left-hand boundary are fixed and the in- 
flow of the plasma is initialized at a constant value. In order to 
trace the condensation of the streaming material onto the cloud 
a new quantity, "colour", is introduced that is set in each cell to 
the density fraction of hot ISM. At the beginning only the cells 
around the cloud possess a non-zero "colour" of value unity. 
During the calculation this quantity is advected like the others, 
such as mass density or energy density. The Poisson equation for 
self-gravity was solved at each tenth timestep because significant 
density changes of the cloud structure happen on a much longer 
timescale than the dynamical one. The energy equation includes 
heating, cooling and heat conduction: 



de 

at 



+ V 



-PV 



T-A- Vq 



(9) 



Here e denotes the energy density, v the velocity, P the pressure, 
r the heating rate, A the cooling rate and q the heat flux. The 
equation of state for an ideal gas is assumed to be valid: 



P=( 7 -l)e 



with 



7 = 5/3 



(10) 



The cooling function used assumes collisional ionisation equi- 
librium and is a combination of the function introduced by 
Bohringer & Hensler (1989) for T > 10 4 K and solar metal- 
licity and by Dalgarno & McCray ( 119721 1 for the lower tempera- 
ture regime. The heating function considers cosmic rays (Black 
1987), X-rays and the photoelectric effect on dust grains (de 
Jong [T977l : de Jong et al. [19801) . 

The heat flux is calculated by taking both the classical and the 
saturated flux into account. In order to apply a smooth transition 
between the classical and saturated regimes we use the analytical 
form by Slavin & Cox ( fl992l ) 



<?sat| 1 ~ exp 



9class 



lQsat| 



(ID 



This guarantees that the lower flux is taken if both differ sig- 
nificantly. The heat flux due to electron diffusion is calculated 
separately using an implicit method which follows a scheme in- 
troduced by Crank & Nicolson (1947) and Juncosa & Young 
( 119711 ). To couple the two directions in space the method of frac- 
tional steps by Yanenko ( 119711 ) is used. 

A detailed description of the implementation of heat conduction 
in the existing hydro-code is given in [Paper I] as well as a com- 
parison of analytical solutions and numerical test cases to prove 
the reliability of the code. 



The initial temperature and density profiles of the clouds are gen- 
erated for hydrostatic and thermal equilibrium under the con- 
straint of spherical symmetry: 



1 



VP 



T(r) = A(r) ; 



(12) 
(13) 



$(r) is the gravitational potential. Setting the temperature Tjsm 
and particle density njsM of the hot and tenuous outer medium, 
the energy density eisM of the plasma is given. The density and 
temperature profile of the cloud is then calculated by integrating 
equations ([T2l and (Qjj) from inside-out using the core temper- 
ature of the cloud as a boundary condition and truncating the 
cloud's outer border where the energy density reaches ejsjyi. 

3. Models 

Here we present the evolution of three different models. Their 
parameters are given in Table [2] For all models the temperature 
of the HIM is fixed to 5.6 • 10 6 K. In two models their density 
is set to 6.6 • 10 -4 cm~ 3 to allow for comparison with static 
models of Paper I and for one model (E) is increased by one 
order of magnitude. 

In model U we consider a massive cloud of 6.4 • 10 4 Mq repre- 
senting a GMC or PCC. This cloud serves as a reference model 
for Paper I in which the fate of this cloud was investigated in an 
identical but static environment. Model E with a cloud mass of 
486.7 Mq and a cloud radius one order of magnitude smaller 
than model U represents a small molecular cloud. Clouds of 
this type can be found as remnants of larger clouds in galac- 
tic chimneys such as the one associated with the Hll region W4 



(Heyer et al. 1996 Taylor et al. 1999[) or behind shock fronts of 



supernovae. The size of model K is similar to model E but the 
mass is decreased to 15.9 Mq to reach an almost homogeneous 
density distribution that can be compared in its evaporation rate 
with Eq. (6) which is valid for uniform clouds. This represents an 
extreme case of a small cloud that is only slightly gravitationally 
bound. 

All models show a typical multi-phase structure. While Model 
U and E possess a dense core and a density decrease outwards, 
Model K is homogeneous. Their radial density profiles are plot- 
ted in Fig.Q] With respect to dealing with more realistic interstel- 
lar clouds with self-gravity and non-equilibrium boundary con- 
ditions, these models differ clearly from the situation implied by 
CM77. 

For comparison all models were studied with and without heat 
conduction. The parameters uo, Ap and T eva are therefore only 
defined for calculations with heat conduction. In all models the 
hot gas streams subsonically with vtsm = 0-3 Mach and all 
clouds are kept at rest so that the relative velocity v le i impor- 
tant for the KH instability (see Appendix B) is identical to vtsm- 
Such a large velocity difference between clouds and the HIM 
is only observed for HVCs, PCCs, and interstellar clouds over- 
taken by supernova shocks. Here we wish to study dynamical 
effects on clouds in addition to heat conduction. For smaller rel- 
ative velocities, as GMCs move through the ISM, the models 
will approach those of numerical heat conduction for clouds at 
rest in a static HIM (presented by us in Paper I). 
The cloud masses were chosen to cover a range between 10 1 
- 10 5 Mq. Because the initial model has to fulfill Eqs. ( fl~2l ) and 
(TTJl ). mass and radius are fixed by the initial central cloud density 
and external pressure and temperature. 
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Table 2. Model parameters used in simulations as described in the text. Valid for all simulations: Tism = 5.6 • 10 6 K i.e. sound 
velocity cjsm = 357 km s -1 , -utsm = w r el = 0.3ctsm = 107 km s _1 



Model 


"ISM A™ 3 


7? cl / pc 


density 
contrast 


M lowr T / M 
(Ea.lB.12l 


M cl /M 


Mmax/ Mq 
(Eq.CJ 


O"0 

(Eq.0 


(Eq.g 


U 


6.6 • itr 4 


41.3 


1.2 ■ 10 4 


2.03 • 10 4 


6.4 ■ 10 4 


1.12 ■ 10 5 


10.5 


1.4 ■ 10 3 


E 


6.6 ■ itr 3 


2.19 


6.1 ■ 10 4 


278.7 


486.7 


1.3 ■ 10 4 


19.8 


2.6 ■ 10 3 


K 


6.6 ■ 10~ 4 


4.0 


5.9 ■ 10 3 


3.08 • 10 5 


15.9 


8.3 • 10 s 


108.4 


1.4 ■ 10 4 



I 



-20 



-,-24 



I 



1 



-26 



-28 



Model U 
Model E 
Model K 



Table 3. Timescales for the three models at their initial configu- 
ration 



30 



40 



Fig. 1. Density distribution for the various initial models. Model 
U and E show a typical multi-phase structure whereas model K 
has a uniform density distribution. 



The saturation parameter <tq (Table 2) indicates a moderately 
saturated heat conduction for model U and E while model K lies 
nearly in the suprathermal regime (Balbus & McKee, 1982 ). The 
importance of heat conduction in these simulations is also indi- 
cated by the fact that the Field length Ap is much larger than the 
cloud radius R c i which means that evaporation of the clouds is 
expected and that the temperature structure of the gas is domi- 
nated by conduction (Begelman & McKee 1990). 
The different dynamical timescales defined as the sound 
travel time over one cloud radius, the timescales t^h for the 
growth of KH instability (see Appendix B) and for evaporation 
are listed in Table|3] Because of the huge density contrast and the 
large gravitational acceleration of the models U and E, KH insta- 
bility should be suppressed, at least at the beginning. Because of 
the stability reasons mentioned in Appendix B and because the 
cloud mass of model E is close to the lower stability limit, this 
cloud is expected to be disrupted during the calculation. Model 
U is further inside the mass limit than Model E and will there- 
fore survive longer than the less massive cloud. Model K on the 
other hand should develop KH instability after tkh = 0.45 Myr. 
Although the evaporation time of Model U and E is too large to 
be reached during the simulation, we should be able to follow 
the disintegration process due to evaporation of Model K. 
In all calculations of T eva the effect of saturation is considered. 

3. 1. Large, massive cloud (model U) 

The inclusion of heat conduction tends to stabilize massive and 
large clouds as illustrated by comparing the density contours at 
the same times (25 Myr, 50 Myr and 75 Myr after the beginning 
of the calculations) for the model U without (Fig. fSJl and with 
heat conduction (Fig. [3). Without heat conduction the edge of 





r dyn = R d/ C c\ 


r KH 


Teva = M/m 


Model 


Myr 


Myr 


Myr 


U 


16.2 




279 


E 


1.95 




77.1 


K 


0.84 


0.45 


3.56 



the cloud is torn by KH instability shortly after the beginning 
as mentioned in Appendix B. A very complex velocity struc- 
ture with many vortices behind the cloud forms. The increase in 
the maximum velocity during the calculation is due to accelera- 
tion of the gas by Bernoulli's effect when it streams around the 
formed cloudlets. It is therefore only a local effect which does 
not affect the large cloud as a whole. In the presence of cool- 
ing stripped-off cloud material becomes thermally unstable and 
forms cloudlets that hit the cloud from the rear. As mentioned by 
Dinge ( 119971 ) these cloudlets produce soundwaves which prop- 
agate through the cloud. On the front side these waves produce 
distortions that can act as seeds for RT instability that lead to 
an elongation of the cloud perpendicular to the stream direction. 
This effect is artifically intensified by the applied numerical sym- 
metry. 

Material once collected on the z-axis sticks there and cannot be 
removed because the velocity in the radial direction is set to zero 
at the axis. A detailed analysis of this artefact and the influence 
of this axis effect is given in Appendix A. 
The cloud develops a compound envelope, with a radially de- 
creasing density in the core region and a diluted outer part of ma- 
terial that is only slightly gravitationally bound. The loss of gas 
by stripping and the recapture from the rear leads to a complex 
mass-loss function with unsteady losses and gains. The mass of 
individual cloudlets that escape into the ISM lies in the range of 
some ten solar masses. 

Fig. 2h shows the evolution of the gravitationally bound cloud 
mass of the simulations without and with heat conduction. For 
comparison the analytical result of the cloud mass is shown tak- 
ing mass-loss by evaporation into account (Eq. 6). From this fig- 
ure we see: 1) The pure evaporation effect on a static homoge- 
neous cloud from CM77 is the most destructive because of its 
shortest timescale, but it overestimates the strength of heat con- 
duction and its heating effect. 2) The pure dynamical destruction 
by KH instability happens on a much longer timescale. As in a 
static self-gravitating and cooling cloud model (Paper I), satu- 
rated heat conduction yields stabilization against evaporation; 
this more realistic approach also reduces the dynamical mass 
loss. While in the static models this is solely caused by an in- 
crease of heat capacity due to energy transfer to denser regions, 
the dynamical destruction is reduced by weakening the KH ef- 
fect. This also means that a cloud that is stabilized in the static 
case by heat conduction is ablated solely and less effectively by 
KH instability. 
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Fig. 2. Evolution of the density distribution for model U without heat conduction. The evolution is shown at the times 25 Myr 
(upper panel), 50 Myr (middle) and 75 Myr (lower). Arrows indicate gas velocities scaling linearly with respect to the maximum 
velocity shown in the upper left. The greyscale represents the density distribution (logarithmic scale). The contour lines represent 5, 
10, 50, 100,. . . x pism- The thick contour encloses the gravitationally bound part of the cloud. 



Obviously the dynamical interaction between the gas phases is 
the dominating effect for cloud destruction in all models. 

After 81 Myr the cloud mass with heat conduction reaches a 
value approximately 7% larger than without, indicating that heat 
conduction reduces the ablation of material from the cloud. A 
comparison of Fig. [2] with Fig. [3] shows significant differences 
in both the cloud structure and its evolution. Heat conduction 
suppresses large-scale KH instabilities as indicated by the al- 
most laminar flow pattern around the cloud. During the whole 



simulation only a single large circulation in the slipstream of 
the cloud is visible. This smoothing process can be understood 
by conducting electrons moving along the steepest tempera- 
ture gradient, i.e. radially toward the cloud, so that the stream 
close to the contact interface is decelerated, the velocity gradi- 
ent decreases and the transition layer (described in Appendix B) 
grows. Numerically, the transition layer is by definition the re- 
gion in which the velocities transit from the HIM to the clouds' 
thermal value. 
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Fig. I3J? shows that over the entire simulation time the thickness 
of the transition layer perpendicular to the stream is, in the heat 
conduction case, nearly twice as thick as in the non-conductive 
case. Comparison with Fig. |B.2| reveals that the transition layer 
without heat conduction lies in the instability zone while it ap- 
proaches or even exceeds the critical layer thickness d for sta- 
bility. This effect is also responsible for the more efficient ac- 
cretion and incorporation of streaming material into the cloud 
visible in Fig. HJ;. The lower velocities in the contact interface 
enable the streaming hot ISM to condense onto the cloud sur- 
face. Thereafter it is very efficiently mixed with cloud material 
and transported into deeper layers of the cloud by turbulent mo- 



tion. This turbulence is also driven by heat conduction that acts 
as an additional energy source in low-density regions where it 
exceeds the cooling. Without heat conduction the accreted ma- 
terial is peeled off from the surface by the grazing high-velocity 
gas before it could be incorporated into the cloud. 

The multi-phase cloud in our simulation can be divided into two 
parts (see Fig. 5): 1) regions smaller than Ap where the heat 
conduction exceeds the cooling so that the density distribution 
becomes homogenized; 2) very dense regions, such as near the 
core where the density structure is unaffected by heat conduc- 
tion. A further surprising result from the numerical models is 
that the heat conduction becomes less effective after the low- 
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evaporation into account, (b) Thickness of the transition layer perpendicular to the stream. The thick solid and dashed lines represent 
the mean value for both cases, (c) Evolution of the amount of accreted material, (d) Evolution of the mean density. 



density outermost shells of the cloud are dissolved. The resis- 
tance of the major part of the cloud against the disruptive KH 
instability and evaporation processes can be revealed from the 
temporal behaviour of the mean density of the bound mass (see 
Fig. Si). Until 2. 

S^dyn m e decline of the mean density is due 
to the dissolution of the outer parts of the cloud. Thereafter the 
mean density increases again and remains nearly constant at a 
value of around 80% of the initial mean density. In contrast the 
decrease of the mean density in the model without heat conduc- 
tion shows the continuous dispersion of the whole cloud by low- 
ering the density also in the core regions. This leads to a flat 
gravitational potential which further destabilizes the cloud. 



3.2. Dense cloud (model E) 

To investigate the effect of heat conduction on denser, less mas- 
sive, i.e. more compact clouds we increase the external density 
PlSM by one order of magnitude. Because the evolution of this 
model is not as violent as the previous one, a snapshot of the evo- 
lution 22.5 Myr after the begining of the calculations is shown 
in Fig. [5] without (upper part) and with heat conduction (lower 
part) for comparison. 



In the non-conductive case even after 22.5 Myr no KH instabil- 
ity is visible because of the stabilizing effect by the high gravi- 
tational acceleration at the cloud surface. 

While the shape of the cloud remains unaltered, variations of the 
mean cloud density and the specific internal energy are present in 
Fig. |6t. These pulsations stem from quasiadiabatic oscillations 
with a timescale given by r a( j = -^3^/ [(37 — 4)Gp\. p is the 
mean cloud density. 

For model E and the adiabatic index 7 = 5/3 , r at i is 7 Myr 
which agrees well with the obtained one. These oscillations orig- 
inate from the interaction of the plasma stream with the cloud at 
the beginning of the simulation. This process triggers accous- 
tic waves that run through the cloud that lead locally to density 
enhancements. 

In addition, the static models ( |Paper IJ > demonstrate that heat 
conduction produces inward travelling sound waves which are 
able to transport energy from the conductive interface to the 
denser cloud interior. By this additional effect the cloud releases 
the transferred energy. 

In the case with heat conduction, small-scale mixing at the cloud 
surface with the outer HIM leads to an increase of the cloud vol- 
ume visible as a decrease in the mean density (Fig. ®>). The 
oscillations occurring in the non-conductive case are efficiently 
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damped by turbulent motions of the cloud material driven by 
heat conduction. The incorporation of hot material into the cloud 
acts as an additional heating agent that raises the temperature and 
therefore the pressure of the cloud and results in an expansion. 
This temperature increase is reflected by the increase of the spe- 
cific internal energy (Fig. [6}:). Regions of the cloud near the vor- 
tex in the slipstream are especially heated and diluted because of 
the more efficient mixing there. This material is only slightly 
gravitationally bound and therefore preferentially stripped off 
the cloud which can be revealed by the distorted cloud shape in 
those regions. Also the small mass-loss of 2% at about 25 Myr 
(12.8 T(jy n ) is caused by this effect. The distortion of the cloud 
to a more streamlined shape as visible in the lower part of Fig. [5] 
leads the flow around the cloud in a more laminar manner. 

As in model U the amount of accreted material is again higher 
with heat conduction than without (Fig. [6}i). This is caused by 
the suppression of the velocity differences in the transition layer. 
However, in the conductive model the absolute accretion rate 
normalized to the total mass is much higher here on absolute 
timescales (after 20 Myr e.g. 0.3% are accreted in model E at 
constant rate but almost 0.1% in model U) but lower on Tdy n ( 
0.26% within 10 Td yn ,E m model E, while 0.25 % within 3 T(jy n jj 
in model U). 



In order to understand this effect quantitatively for the different 
parameters of the two models one can make a zeroth-order ap- 
proach neglecting density and mass changes. An accretion rate 
il can be defined as Af acc per dynamical time Td yn and reformu- 
lated by normalization to the total mass Mtot as shown for model 
U in FigHJ; and for model E in Fig|6jl, respectively, by 



n = 



Ma, 



Mtr 



Mtot T, 



dyn 



(14) 



where Td yn must not be further determined. Applying this rela- 
tion, the ratio 57tj / fig follows as 

n v = 0.0025 M cUJ 0.0026 M cl;E ) ^ 

3 T dyn jj 10 T dyn E 

With Eq. (O, Mace = « r dyn follows so that one can write 
the ratio M acc \j /M acc g where the different Td yn cancel out with 
Eq. ( fTBT l and one obtains the value 421. This value agrees as 
an order-of-magnitude effect within the uncertainties and being 
aware of the simple estimates of the structure, if one assumes 
that the condensation A/ a cc depends on the cloud's surface, the 
surrounding density and the Field length is given by: 



M a , 



oc A F Rl x n ISM — — — n IS M 



(16) 
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This ratio for models U and E is 361. 

Compared to model U where heat conduction prevents the dis- 
ruption of the cloud, in this model the cloud is already stabilized 
without heat conduction and condensation dominates the evo- 
lution. The dip in the mass loss curve (Fig. 9a) after 18 Myr 
is caused by the separation of a small cloudlet with successive 
rebounce so that it remains invisible in Fig. 9d. Evaporation, 
however, is then not directly caused by heat conduction. Instead, 
conduction expands the cloud and KH instabilities dissolve the 
back parts of the cloud. Fig. [6^ reveals that the total bound 
mass remains constant for the case without heat conduction but 
decreases with heat conduction only by 2% within 12.8 r^ yn . 
Interestingly both conductive models U and E experience an ex- 
pansion that reduces the mean density by nearly 20 %. Although 
heat conduction leads to a slight expansion of the cloud and, 
by this, to this minor fraction of mass loss, the analytical result 
in contrast requires evaporation of 30% of the cloud after 12.8 
r dyn- The analytical prescription overestimates the real effect by 
far, when both gas phases are dynamically different. 

3.3. Small, homogeneous cloud (model K) 

A third model describes an extreme case: Model K represents a 
small cloud with no distinct core and a density contrast of almost 



one order of magnitude lower than the stable model E. The an- 
alytical consideration of the classical heat conduction (CM77), 
therefore, requires fast and efficient evaporation (see Tableland 
Fig. [9J3). This cloud is only marginally gravitationally bound, 
so that its gravitational acceleration cannot sufficiently stabilize 
it against KH instability (see Appendix B). The timescale for 
the onset of KH instability is approximately tkh = 0.45 Myr, 
i.e. almost half of id yn . The cloud should therefore be destroyed 
within a few dynamical timescales. The violent action of KH in- 
stability is easily visible in Fig. [7] where the density contours for 
Model K without heat conduction are shown in steps of 3.4 Myr 
after the beginning of the calculations. The density contours in 
the case with heat conduction are revealed in Fig.[8]for the same 
timesteps. 

Fig. [7] shows that the cloud model K is already very fragile at 
the beginning of the simulation. The cloud is soon perpendic- 
ularly stretched because of the low pressure at both sides due 
to Bernoulli's theorem. The density in these regions is lowered 
and its cooling ability diminished. This, together with the mix- 
ing of hot ISM into the cloud due to small scale perturbations 
triggered by means of KH instability leads to an increase of the 
temperature which can be discerned as an increase of the spe- 
cific internal energy (Fig. [9}?). This energy input is sufficient to 
unbind the cloud material. From Fig. [9]j one also sees that the 
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Fig. 7. Evolution of the density distribution for the simulation K without heat conduction. The evolution is shown at the times 
3.4 Myr (upper panel), 6.8 Myr (middle) and 10.2 Myr (lower). Arrows indicate gas velocities scaling linearly with respect to the 
maximum velocity shown in the upper left. The greyscale represents the density distribution (logarithmic scale). The contour lines 
represent 5, 10, 50, 100,. . . x pism- The thick contour encloses the the domain where p c \ > 100 x pism- 



stripped-off gas contains a larger specific internal energy than 
the rest, i.e. energy that is acquired from dynamically excited 
turbulence is lost. Due to the initiated turbulent velocity field in- 
side the cloud the mass distribution changes throughout the cal- 
culation. Consequently, no distinct core can be formed. The low 
binding energy allows the cloud material to be caught by the gas 
flow and to be carried away at about 5 Myr (5.9 T^y n ). Because 
the stripped cloud material is disrupted into small packages, the 



mass-loss rate (Fig. [9^) shows a step-wise function. After about 
10.8 Myr (12.8 Tdyn) the cloud is nearly destroyed and only 30% 
of the initial cloud mass has survived. 

The evolution in the case with heat conduction is very similar 
during the first 5 Myr to that without heat conduction. Both the 
elongation of the cloud and the mass-loss rate show nearly the 
same behaviour. In contrast to the non-conductive case, KH in- 
stability is suppressed and thus the additional heat input through 
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Fig. 8. Same as Fig. [6] Model K with heat conduction. 



mixing with the hot ISM is reduced. The temperature and there- 
fore the specific internal energy rises only moderately (Fig. [9}}). 
Consequently, the velocity field within the cloud is less turbulent 
so that a core can form. The cloud interior develops a relatively 
stable mass distribution with a dense core which stabilizes this 
model against mass loss at least twice as long (Fig. |9^). Due 
to heat conduction a transition layer with a radial decrease of 
the density and velocity gradient has formed and with this a ra- 
dial temperature gradient along which the cloud material flows. 
Some material is lost after 7.5 Myr when a large cloudlet with a 



mass of about 30% of the whole cloud decouples from the elon- 
gated cloud at its top. 

As in large massive clouds (model U) where the development of 
a transition zone near the cloud surface leads to a reduction of 
KH instability, the same is valid for small homogeneous clouds. 
No large-scale KH instability is seen although the initial state of 
the cloud is very unstable against it. At the end of the calculation 
the total mass of the cloud is about 30% higher in the case with 
heat conduction than without. 

The heat conduction process leads to dynamical processes at the 
clouds' surface that decelerates the disintegration of the whole 
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cloud by suppressing the large scale KH instability. In addi- 
tion, considering cooling in the heat conduction description sta- 
bilizes this cloud against mass loss, while the analytical solu- 
tion would again destroy the cloud by evaporation much more 
rapidly, within 6 T^ yn . 

4. Discussion 

In agreement with analytical results for the static case (CM77) 
we have proven that heat conduction must not be neglected in 
investigations of the evolution of interstellar clouds in a hot di- 
lute plasma and, in general, of the coexistence of the multi -phase 
ISM. While CM77 indicate that all the clouds in this paper would 
undergo evaporation due to heat conduction, we demonstrated in 
|Paper I| and in sect. 1 .3 that a realistic approach that accounts for 
the electron conductivity properly leads to the opposite results, 
namely, that condensation dominates. The cloud conditions ap- 
plied to our models differ from the consideration of CM77. In 
this paper we have shown with models that the fate of a cloud 
and its evolution is drastically changed by a gas stream and de- 
pend strongly on the mass and its internal structure of the cloud, 
i.e. mainly on its binding energy. In general, the relative motion 
of a subsonically hot plasma stream stabilizes the clouds. 
Without heat conduction, clouds with their initial states close to 
or inside the KH-unstable regime suffer from huge mass loss in 
the form of stripped-off cloudlets. While a small homogeneous 
cloud (model K) is stable for about 5 T^y n and then strongly ex- 
posed to disruption into small gas packets, small dense clouds 
(model E) can avoid the transition into the KH-unstable state 
and resist the violent hot plasma so that mass loss or even any 
strong deformation of the cloud does not occur. Large massive 
clouds (model U) lose about 17% of their mass within 5 T& yn and 
may dissolve on larger timescales. 

For the static case of the coexistence of two gas phases, namely, 
clouds embedded in a HIM, the mass transfer according to heat 
conduction can be calculated analytically (CM77). We have cho- 
sen plasma conditions that would require evaporation of cloud 
material and therefore mass loss from the cloud from analytical 
considerationsly. Under the dynamical action of a relative mo- 
tion between the gas phases the conditions change with respect 
to the static case in two ways: Dynamical instabilities can change 
the shape of the cloud and can increase its surface. In contrast to 



the analytical results, the state of the hot ISM remains constant 
because of its replenishment by the fixed streaming conditions 
and therefore cannot react to the evaporation and condensation 
process and by this e.g. self-regulate the mass transfer to find an 
equilibrium (see e.g. Koppen et al. 1998). 
Large and massive clouds survive longer in the hot plasma flow 
with heat conduction than without. Because of electron invasion 
through the surface, a transition zone forms at the edge of the 
cloud where density and velocity gradients are lowered. A state 
can be reached where the KH instability is suppressed and the 
formerly unstable cloud becomes stabilized. This can be shown 
analytically and numerically reproduced (model U). Since the 
evaporation rate is much less than the one predicted by ICM77I 
these facts lead to a cloud mass at the end of the calculation 
that is even slightly higher (7%) than in the case without heat 
conduction. 

Although the maximum cloud mass implied here is much 
lower than those of HVC complexes moving through the galac- 
tic halo with masses of a few 10 6 Mq such as Complex C 
(Wakker et al. 1988) or compact HVCs located in the intergalac- 
tic medium dBraun & Burton 1999K only clumpy substructures 
seem to decouple from complexes and approach the galactic 
disk and experience on their path through the halo interaction 
with the hot gas that leads to the observed head-tail structures 
(Brims et al. 2000). Their masses, on the other hand, range from 
a few umpteen solar masses to a few 10 4 M ^, like recently 
found compact HVCs in the inner galaxy dStil et aL~,2 006). 
Nevertheless, calculations with even higher masses are in prepa- 
ration but it can be expected that the tendency to stabilize the 
cloud and to reduce the ablation of material from the cloud will 
continue. Heat conduction is therefore a physical process that 
enhances the dynamical stabilization and has to be taken into 
account in the consideration of the survival of HVCs. 
For the PCCs, heat conduction offers a mechanism to incor- 
porate metal-rich hot gas that becomes homogeneously mixed. 
Even with a low accretion fraction of only 0.3% a hot gas 
metal content of solar and above, which is reasonable e.g. 
from X-ray determinations of the halo gas around giant ellip- 
ticals dMatsushita et al. 20031 , would lead to almost 1/100 Zq. 
Globular clusters formed from PCCs and enriched by this accre- 
tion mechanism caused by heat conduction must be expected to 
show a large range of metallicities. 
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When star formation sets in, all protostars are formed from 
molecular clouds with nearly equal metallicity. The absence of a 
significant spread in [Fe/H] in most globular clusters (Freeman 
& Norris flMTl Fahlmann et al. [T985l Norris [198"8l is an indica- 
tion that the stars have formed out of well mixed metal-enriched 
substrates that could have been polluted by an external source 
(Murray & Lin 1 1990t . This mechanism provides an alternative 
explanation to the self-enrichment scenario of globular clusters 
dBrownetal. 19911 [19951 . 

When looking at smaller clouds one has to distinguish between a 
homogeneous (model K) and a clearly centrally peaked density 
distribution (model E). In the latter case the gravitational poten- 
tial is strong enough to stabilize the cloud against large-scale 
perturbations triggered by the dynamical action of the streaming 
ISM. The influence of the HIM is limited to an additional heat 
input due to heat conduction and small-scale mixing especially 
in regions near the vortex in the slipstream of the cloud. Because 
of the high density in the core regions, the additional heat input 
is nearly compensated by cooling. Only the rim of the cloud be- 
comes diluted and so susceptible to evaporation. Thus, the mass 
loss (2%) is reduced by more than a factor of 15 with respect 
to the analytical approach (~ 30%). In this case the fate of the 
cloud depends only on small-scale, not on large-scale dynamics. 

These results cannot be extrapolated to homogeneous clouds 
(model K), where already a small additional heat input is suf- 
ficient to unbind the whole cloud. Although KH instability is 
also suppressed in this case, the cloud becomes elongated by the 
Bernoulli effect, which is very efficient because of the flat gravi- 
tational potential. It is this large-scale phenomenon that peels off 
a large fragment of the cloud. Evaporation due to heat conduc- 
tion is negligible in comparison to the mass loss due to the large- 
scale effects. Heat conduction acts instead as an agent to extend 
the survival time of the cloud by the suppression of KH instabil- 
ity and therefore the reduction of the mass-loss in comparison to 
the non-conductive case. While the analytical approach of model 
K leads to cloud destruction by evaporation within almost 6 T&y B 
(Fig. [9^), in the numerical simulation the dynamics stabilize the 
cloud and heat conduction extends the destruction time to more 
than 1 1 T(iy n . After the loss of a large fragment with 30% of the 
cloud mass, the total disruption of the cloud is most likely a con- 
sequence of the dynamics of the streaming ISM. It is therefore 
not surprising that the analytical evaporation rates of CM77 are 
incompatible with our calculations. 

Comparing the CM77 evaporation rates with computed models 
for the static case (Paper I) reveal that condensation may occur 
for large clouds in temperature regimes where lCM77l also predict 
evaporation. 

The three presented models can only give a first insight into the 
fundamental importance of dynamics and additional heat con- 
duction in the evolution of the ISM phases. Further investiga- 
tions are necessary to determine the dependence of evaporation 
and condensation on the physical state of the phases. This is nec- 
essary for the simulation of galaxy evolution (see e.g. Samland 
etal. 1997, Hensler 1999). 

The effect of magnetic fields on heat conduction must be dis- 
cussed, especially when it is neglected, as in our models. When 
electrons of the HIM enter an interstellar cloud they transfer their 
energy to the cloud by collisions with mainly neutral HI atoms 
of density nni and collisional cross section qni- Their collisional 
mean free path A Cj Hi/cm = (nnr qm) -1 is about 10 16 cm. On the 
other hand, magnetic fields force charged particles with mass 
m e moving with velocity u e jcM perpendicular to the B field 



vectors due to the Lorentz force to gyrate with the so-called 
(Larmor) radius a£ e . For electrons this reads as 



C m e ■ U e jCM 

lei B 



(17) 



Since we assume in our modelling that at the clouds' surface, 
both gas phases achieve pressure equilibrium we can also set 
the magnetic pressure P mag = B 2 /8ir in equilibrium. Replacing 
m e ■ ICM by the ICM temperature Tjcm> it cancels with that 

of the pressure and one gets a^ ^ = 4.46 • 10 5 n Ii 1 J // ( ^ M . 
The ratio of mean free path to electron Larmor radius then is 



A 



= 3.5 • 10" 



c.HI 



i 1/2 



(18) 



For weaker B fields than in gas pressure equilibrium the ratio can 
become much larger, i.e. the free electron motion is less ham- 
pered. But this only holds for the fraction of electrons moving 
perpendicular to the B vector, while electrons moving parallel to 



the B field are unaffected. Malyshkin & Kulsrud (2001) how 



ever, showed that random magnetic fields also reduce the diffu- 
sivity of electrons traveling along the B field lines. This reduc- 
tion of free electron motion and thus heat conduction by mag- 
netic fields can also be parametrized by $ s as in eq|5] 

Appendix A: Axis effects (model U ) 

The numerical models in this paper are calculated assuming 
cylindrical symmetry. This method has numerical artefacts near 
the symmetry axis. Material, once fallen onto this axis, cannot be 
removed from it because of mirrored conditions, i.e. zero gradi- 
ents. As a result a nozzle is formed at the front of the cloud, 
which Dinge (1997) interprets as RT instability. On the other 
hand, ablated material that is pushed towards the axis by vor- 
tices behind the cloud builds up elongated spurs. In contrast, in 
3D simulations these clumps of a typical thickness of a few cells 
oscillate in the slipstream of the cloud because of asymmetrical 
instabilities. In order to test the influence of these artefacts the 
evolution of the clouds we switched to a 2D Cartesian grid and 
so removed the symmetry axis. Fully 3D simulations with suf- 
ficient spatial resolution including heat conduction are still too 
time consuming to be produced. The initial model of the sim- 
ulated clouds are therefore no longer spheres but infinite cylin- 
ders. As a consequence, the velocity field, the gravitational po- 
tential and the density distribution of the initial model differ from 
the axisymmetric case but are close to model U. 
The evolution of this model U* is shown 25 Myr, 50 Myr and 75 
Myr after the beginning of the calculations with consideration 
of heat conduction (Fig. IA.lt . Because of the infinite cylinder 
instead a sphere, the velocity of the stream near the top and the 
bottom of the cloud is much higher than in the simulations shown 
before. As a result, the onset of KH instability is facilitated. The 
tail of ablated cloud material at the rear of the cloud can wave 
in the stream. Regions of the tail that become too dilute are cut. 
Nevertheless a massive tail similar to model U forms. No spur 
forms at the front of the cloud and must therefore be identified 
as an artefact of the axisymmetric conditions. 
Although there are some differences in the appearance of the 
clouds in cylindrical and Cartesian geometry, the global forma- 
tion of a head-tail structure of the cloud remains the same in both 
geometries. The thin elongated structures stripped off the cloud 
in the Cartesian grid are comparable to the material that get stuck 
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Fig. A.l. Evolution of the density distribution for the 2d simulation using a cartesian grid with heat conduction. The evolution is 
shown at the times 25 Myr (upper panel), 50 Myr (middle) and 75 Myr (lower). Arrows indicate gas velocities scaling linearly with 
respect to the maximum velocity shown in the upper left. The contour lines represent 5, 10, 50, 100,. . . x pism- The thick contour 
encloses the gravitationally bound part of the cloud. 



on the symmetry axis in cylindrical symmetry. In both cases the 
filaments leave the computational domain and do not influence 
the evolution of the gravitationally bound cloud. Therefore the 
use of a cylindrical symmetry is justified. 

Appendix B: Stability analysis 

B. 1. Rayleigh-Taylor instability 

If a cool, dense cloud moves through a hot, tenuous gas, its sur- 
face will become subject to KH and RT instabilities. If the cloud 



is massive enough so that the effective gravity is directed toward 
the cloud center, RT instabilities can be suppressed. The growth 
rate for the RT instability in the absence of a gravitational field 
is expressed as (Chandrasekhar [196 1 I t 

I 2 1 ; Pel - PlSM IU 

\uj\ RJ = ak ■ , (B.l) 

Pel + PlSM 

where k is the wavenumber of the perturbation, a the accelera- 
tion of the cloud due to the wind, p Q \ and pism are the densities 
of the cloud and of the streaming medium, respectively. This 
means that perturbations with shorter wavelengths grow faster 
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than larger ones. As long as the wavelength is small compared 
to the cloud radius, the net effect of RT instability is a small- 
scale mixing of the ISM with the cloud but does not lead to its 
disruption. The cloud is stabilized against RT by a gravitational 
field and RT is totally suppressed if the gravitational accelera- 
tion g is larger than the acceleration due to the wind (Murray et 
al. 119951 . 

B.2. Kelvin-Helmholtz instability 

In order to study the influence of KH instability we perturb the 
incompressible fluid equations 



dvi 



-P9i 



dp 
dx 



dp 
dt 

and 

3 

E 



5> 

i=l 



dxi 



= 



[i = 1,2,3] (B.2) 



(B.3) 



dxi 



= 



(B.4) 



for a solenoidal velocity field according to 
Chandrasekhar (1961). The quantities density p, pressure 
p and velocity v, are replaced by the perturbed ones p + 6p, 
p+Sp and v = Vq + 8v. The acceleration due to the gravitational 
force g remains unperturbed. The initial velocity field Vq for the 
KH stability analysis is zero except for v x (z) which is the flow 
velocity and therefore an arbitrary function of the height z. After 
linearization we have five partial differential equations (PDEs) 
for the five perturbed quantities. We now consider perturbations 
of the form 



For unstable wavenumbers the timescale for the exponential in- 
crease of the perturbation reads 



1 



TKH 



(B.9) 



Because most destructive perturbations are those with wave- 
lengths of the order of the cloud radius (A = R c \ = 2n/k) there 
is a critical gravitational acceleration g cr ; t above which the cloud 
is stabilized: 



2™relPlSM 
Pcl-Rcl 



(B.10) 



If a spherical cloud with uniform density is assumed, its gravita- 
tional acceleration is 



GM, 



cl 



and M d = 4/37rp d i^ 



This leads to a critical mass of the cloud 



Mio 



V6n Vl 



( Pism\ 
Pel \GpciJ 



3/2 



(B.ll) 



(B.12) 



above which KH instabilities are suppressed. Clouds with 
A/i ower < M c \ < Af max should survive long enough so that star 
formation can occur. These estimates are valid for the case of 
density and velocity jumps at the cloud boundary. 
If we take dissipative forces like heat conduction or viscosity 
into account, which are able to diminish the density and velocity 
gradient, the picture changes drastically. For further investiga- 
tions we introduce an artificial transition layer of thickness 2d 
and density po in which the velocity decreases linearly from the 
value Uo outside the cloud (z > d) to — Uo inside the cloud 
(z < —d) (see Fig. IB.lt . The ISM at velocity Uq has a density 
po ■ (1 — e), the cloud material with — Uq has po ■ (1 + e). 



SX oc exp[i(k x x + k y y + uit)] . 



(B.5) 



Inserting this into the five PDEs and combining them we have 
a differential equation for the z-dependence of the perturbation 
6v z which is perpendicular to the streaming v x 



d 

dz 



p(u> + k x v x ) 



dSv, 
dz 



-k 2 p8v z (uj + k x v x ) = gk 2 ( — 



Sv z 



dz / lu + k x v x 



(B.6) 



This equation can be solved analytically for two homogeneous 
fluids of densities pism an d Pel separated by a horizontal bound- 
ary and moving with horizontal velocities wism and v c \. Starting 
with a stable configuration, i.e. the density of the upper fluid 
PlSM i s l ess m an the lower one, we obtain the dispersion relation 
forw(fc) 



d - 

- 
-d 



density velocity 



Po(!-e) 



p (l+e) 



U z/d 



-U U 



HIM 



transition 
layer 



cloud 



Fig. B.l. Density and velocity structure for the stability study: 
a cloud moving with 2Uo relatively to a hot gas with density 
contrast (1 + e)/(l - e). 



I k 2 p c ipi SM (v c i - vism) 2 - fffcQci - Pism) 



(Pel + Pism) 2 



(B.7) 



Instability occurs for real values of u(k), so there is a mini- 
mal wavenumber fc m ; n above which no stable mode exists. With 
u rel = v cl — V ISM an d Pel ^ PlSM this upper limit for stability 
reads: 



3Pcl 



PISM V M 



(B.8) 



For a stability analysis we solve Eq. lB.6l together with the bound- 
ary conditions at the transition zone. The relation between u> and 
k as a function of d and e is given by 



exp[— 2 re] 



I - 



1 - 



re(^ + l) 2 



J+ (v+ 1) + 0.5ere(i/+ l) 2 

K(V- l) 2 



J — [y — 1) -0.5e«(i/- l) 2 



(B.13) 
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with the Richardson number J for this problem 

tad lu 
J=-^T, and n = 2kd. (B.14) 

Modes with wavenumbers k are stable if there are no complex 
solutions for to. The case e — > is discussed in Chandrasekhar 
dl9611 l. For arbitrary e we solved Eq. ( IB. 13b numerically, yield- 
ing an instability strip in the fc-space. For some combinations of 
d and e, most destructive perturbations with wavelengths of the 
order of the cloud radius can be stabilized. 
As an example, we investigate a model with a density contrast 
of 1.2 ■ 10 4 , comparable to the Model U in subsequent calcula- 
tions, which corresponds to e = 0.99983. For different values of 
the transition zone thickness d, the consequence of perturbations 
by certain wave vectors k is calculated and shown in Figure |B~2l 
From this, two facts become obvious: clouds that are unstable 
for certain k are stabilized if the transition zone thickens. In this 
particular case, e.g., perturbations with wavelengths on the order 
of the cloud radius (k = 2ir/R c i) that are usually most destruc- 
tive, are damped for transition layer thicknesses d of more than 4 
pc (upper right zone in Fig. |B.21 >. Secondly, self-gravity leads to 
smaller ds and, by this, to stability for low ks (lower left region 
in Fig. |B.2t . For k = 2tt/R c i this model is stable for d smaller 
than 0.002 pc. 

Spatially poorly resolved numerical simulations produce artifi- 
cially thick boundary layers due to numerical noise. Since such 
extensions would put the model into the instability strip, KH in- 
stability sets in. Heat conduction broadens the boundary layer to 
values d > 4 pc so that the cloud becomes stabilized. 
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Fig. B.2. Kelvin-Helmholtz instability strip for Model U with 
a density contrast e of 1.2 • 10 4 at various wave numbers k and 
extents of the transition layer d. 
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